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Outflows and accretion in a star-disc system with 
stellar magnetosphere and disc dynamo 
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Abstract. The interaction between a protostellar magnetosphere and a surrounding dynamo-active accretion disc is investi- 
gated using an axisymmetric mean-field model. In all models investigated, the dynamo-generated magnetic field in the disc 
arranges itself such that in the corona, the field threading the disc is anti-aligned with the central dipole so that no X-point 
forms. When the magnetospheric field is strong enough (stellar surface field strength around 2 kG or larger), accretion happens 
in a recurrent fashion with periods of around 15 to 30 days, which is somewhat longer than the stellar rotation period of around 
10 days. In the case of a stellar surface field strength of at least a few 100 G, the star is being spun up by the magnetic torque 
exerted on the star. The stellar accretion rates are always reduced by the presence of a magnetosphere which tends to divert a 
much larger fraction of the disc material into the wind. Both, a pressure-driven stellar wind and a disc wind form. In all our 
models with disc dynamo, the disc wind is structured and driven by magneto-centrifugal as well as pressure forces. 
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1. Introduction 

The interaction of a stellar magnetic field with a circumstel- 
lar accretion disc was originally studied in connection with ac- 
cretion discs around neutron stars (Ghosh et al. 1977, Ghosh 
& Lamb 1979a,b), but it was later also applied to protostel- 
lar magnetospheres (Konigl 1991, Cameron & Campbell 1993, 
Shu et al. 1994). Most of the work is based on the assumption 
that the field in the disc is constantly being dragged into the 
inner parts of the disc from large radii. The idea behind this is 
that a magnetized molecular cloud collapses, in which case the 
field in the central star and that in the disc are aligned (Shu et 
al. 1994). This was studied numerically by Hirose et al. (1997) 
and Miller & Stone (1997). In the configurations they consid- 
ered, there is an X-point in the equatorial plane (see left hand 
panel of Fig. [0, which can lead to a strong funnel flow. 

Another alternative has been explored by Lovelace et al. 
(1995) where the magnetic field of the star has been flipped 
and is now anti-parallel with the field in the disc, so that 
the field in the equatorial plane points in the same direction 
and has no X-point. However, a current sheet develops above 
and below the disc plane (see right hand panel of Fig. [Q. 
Numerical simulations of such a field configuration by Hayashi 
et al. (1996) confirm the idea by Lovelace et al. (1995) that 
closed magnetic loops connecting the star and the disc are 
twisted by differential rotation between the star and the disc, 
and then inflate to form open stellar and disc field lines (see 
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also Bardou 1999, Agapitou & Papaloizou 2000). Goodson et 
al. (1997,1999) and Goodson & Winglee (1999) show that for 
sufficiently low resistivity, an accretion process develops that 
is unsteady and proceeds in an oscillatory fashion. The inflat- 
ing magnetosphere expands to larger radii where matter can be 
loaded onto the field lines and be ejected as stellar and disc 
winds. Reconnection of magnetic field lines allows matter to 
flow along them and accrete onto the protostar, in the form of 
a funnel flow (see also Romanova et al. 2002). Consequently, 
their stellar jets show episodic behaviour; see also Matt et al. 
(2002). 

The gas in the disc is turbulent and hence capable of con- 
verting part of the kinetic energy into magnetic energy by dy- 
namo action. Furthermore, the differential rotation of the disc 
can allow for large scale magnetic fields, possibly with spatio- 
temporal order (similar to the solar 11 year cycle; see, e.g., 
Parker 1979). Such a dynamo may have operated in the so- 
lar nebula and in circumstellar or protoplanetary discs (e.g., 
Reyes-Ruiz & Stepinski 1995). However, not much is known 
about the mutual interaction between a stellar magnetic field 
and a dynamo-generated disc field. In particular, we need to un- 
derstand the effect of the magnetic field generated in the disc, 
and we also want to know how this is being affected by the 
presence of the central dipole field. Following a recent attempt 
to model outflows from cool, dynamo-active accretion discs 
(von Rekowski et al. 2003; hereafter referred to as Paper I), 
we present in this paper a study of the interaction between the 
dynamo-generated disc magnetic field and the field of the cen- 
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Fig. 1. Sketch showing the formation of an X-point when 
the disc field is aligned with the dipole (on the left) and the 
formation of current sheets with no X-point if they are anti- 
aligned (on the right). The two current sheets are shown as 
thick lines. In the present paper, the second of the two con- 
figurations emerges in all our models, i.e. with current sheets 
and no X-point. 

tral star. In our model, both the star and the disc have a wind, 
resulting in a structured outflow that is driven by a combina- 
tion of different processes including pressure-driving as well 
as magneto-centrifugal acceleration. 

2. The magnetospheric model 

A detailed description of the model without magnetosphere can 
be found in Sect. 2 of Paper I. We begin the description of the 
present model by reviewing the basic setup of the model in 
Paper I. The implementation of the stellar magnetosphere is 
described in Sect. 1241 



2. 1. Basic equations 

We solve the set of axisymmetric MHD equations in cylindrical 
polar coordinates (zu, ip, z), consisting of the continuity equa- 
tion, 
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the Navier-Stokes equation, 
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and the mean field induction equation, 
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Here, u is the velocity field, g is the gas density, p is the gas 
pressure, $ is the gravitational potential, t is time, D/D< = 
d/dt + u ■ V is the advective derivative, F = J x Z? + V • r is 
the sum of the Lorentz and viscous forces, J = V x £?//io is 
the current density due to the mean magnetic field £?, /io is the 
magnetic permeability, and r is the (isotropic) viscous stress 



tensor. In the disc, we assume a turbulent (Shakura-Sunyaev) 
viscosity, v t ~ assCs^o, where ass is the Shakura-Sunyaev 
coefficient (less than unity), c s = (-/p/ g) 1 ^ 2 is the sound speed, 
7 = Cp/cv is the ratio of the specific heat at constant pressure, 
c p , and the specific heat at constant volume, c v , and zo is the 
disc half-thickness. 

As described in Paper I, in order to maintain a statistically 
steady accretion disc, we need to replenish the mass that is 
accreted through the disc and onto the star. We therefore in- 
clude a mass source, <7g lsc , that is restricted to the disc and 
self-regulatory, i.e. it turns on once the local density in the disc 
drops below the initial density distribution, go(r), of the hy- 
drostatic equilibrium. This means that matter is injected into 
the disc only wherever and whenever g < go, and the strength 
of the mass source is proportional to the gas density deficit. 
Therefore, we do not prescribe the distribution and magnitude 
of the mass source beforehand, but the system adjusts itself. 
We also allow for a self-regulatory mass sink, q^ al , at the po- 
sition of the central object (protostar) to model accretion onto 
the central star without changing the stellar radius. The mass 
sink is modelled in a way analogous to the mass source. In the 
models with magnetosphere, however, a mass sink in the star is 
modelled by setting density, velocity and magnetic field in the 
star to their initial values at each timestep. This is necessary in 
order to anchor the magnetosphere in the star (cf. Sect. 12. 4> . 
The mass source appears also in the Navier-Stokes equation, 
unless matter is injected with the ambient velocity of the gas. 
We always inject matter with Keplerian speed, %. This leads 
to an extra term in the Navier-Stokes equation, (uk — u)qg lsc , 
which vanishes only if the gas rotates already with Keplerian 
speed or when no mass is injected (qg lsc = 0). 

A physically realistic accretion disc is much more strongly 
stratified than what can be represented in the simulations. We 
expect most of the actual accretion to occur in the innermost 
parts of the disc. These inner parts are also much cooler and 
therefore they should be spinning at almost exactly Keplerian 
speed, i.e. usually faster than the outer parts of the disc. This 
is the reason why we choose to inject new matter at Keplerian 
speed. One additional reason why we do not inject matter with 
the ambient gas velocity is that we want to prevent a runaway 
effect (cf. Sect. 2. 1 in Paper I). Such a runaway could result in a 
loss of the initial angular momentum of the disc, which would 
never be replenished, so that the entire disc would eventually 
lose its angular momentum and all matter would be accreted. 
However, in some cases we have compared simulations with 
and without Keplerian injection and found the difference to be 
small. 

Furthermore, mass replenishment in our disc is necessary, 
because we model an accretion disc that is truncated at about 
0.19 AU (cf. Sect. 12.31 . i.e. we model the inner part of a re- 
alistic (protostellar) accretion disc. The total run time of our 
simulations is between ^150 days (Model S) and ~ 1900 days 
(Model Ml), corresponding to about 33 to 422 Keplerian or- 
bital periods of the inner disc edge (cf. Sects 12. 3l andl3l. This is 
shorter than the life time of the inner part of a protostellar ac- 
cretion disc. On the other hand, a typical advection time in our 
models is (ra7 out — n7 in )/u ro w (1.9 — 0.6) x 0.1 AU/ (0.15 x 
10 2 krns" 1 ) w 15 days (cf. Sects I2.3l and l3.5l . which is much 
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shorter than the typical run time so that our disc would dis- 
appear if there were no mass supply. One way to model mass 
accretion from the outer parts of the disc to the inner parts, is to 
extend our disc to the radial boundary of our computational do- 
main and to inject matter at this radial boundary. Another way 
- that we have chosen - is to inject matter locally in the disc 
wherever and whenever it is needed. 

The technique of invoking self-regulatory terms that act 
in certain parts of the computational domain is an alternative 
to prescribing boundary conditions. Similar techniques have 
been used in modelling the Earth's magnetosphere (Janhunen 
& Huuskonen 1993) and in magnetized Couette flow problems 
where inner and outer cylinders are present as part of a carte- 
sian mesh (Dobler et al. 2002). In the present case, our star 
and accretion disc play a hybrid role in that they represent not 
only boundary conditions to the corona including the region be- 
tween disc and star. In the models without magnetosphere, the 
star is not only a region where mass is being absorbed by the 
mass sink, but one can also study the dynamics of the velocity 
and magnetic fields in the star and the resulting effects. We do 
not prescribe inflow at the star's surface; instead, our model al- 
lows for both inflow and outflow at different parts of the star's 
surface. In the models with magnetosphere, the star remaining 
within the mesh allows one to study the effects of a spinning 
dipole. The disc surface is not an outflow boundary condition, 
i.e. a region where mass is being injected into the corona, but 
the outflow develops itself. Resolving the disc makes it possi- 
ble to include a disc dynamo which is the innovative ingredient 
in the models of Paper I and the present paper. Our modelling 
of the disc dynamo will be explained in the following section. 

We assume that the magnetic field in the disc be gener- 
ated by a standard a 2 £l dynamo (e.g., Krause & Radler 1980), 
where a is the mean-field a effect and O is the angular velocity 
of the plasma. This implies an extra electromotive force, aB, 
in the induction equation for the mean magnetic field, B, that 
is restricted to the disc. As usual, the a effect is antisymmetric 
about the midplane with (see Paper I) 

z £disc(r) fA . 

a = a — — — 2~T~2> ( 4 ) 
z 1 + v%/cj 

where is the Alfven speed based on the total magnetic 
field, £disc is a profile specifying the shape of the disc (see 
Sect. 12. 21 . and ao is a parameter that controls the intensity 
of dynamo action. We choose the a effect to be negative in 
the upper half of the disc (i.e. ao < 0), consistent with re- 
sults from simulations of accretion disc turbulence driven by 
the magneto-rotational instability (Brandenburg et al. 1995, 
Ziegler & Riidiger 2000). The resulting magnetic field sym- 
metry is roughly dipolar. This symmetry is seen both in three- 
dimensional simulations of accretion disc dynamos driven by 
turbulence from the magneto-rotational instability and in so- 
lutions of the a£l dynamo problem with the a effect negative 
in the upper disc half, provided that the accretion disc is em- 
bedded in a conducting corona (e.g., Brandenburg et al. 1990, 
Brandenburg 1998, von Rekowski et al. 2000, Bardou et al. 
2001). Further, we include a quenching which leads to the 
disc dynamo saturating at a level close to equipartition between 



magnetic and thermal energies. The (constant) magnetic diffu- 
sivity r\ is finite in the whole domain, and enhanced in the disc 
by turbulence. We do not take into account that the diffusivity 
in the star might also be enhanced due to convection but assume 
the same value in the star as in the corona. 

To ensure that B is solenoidal, we solve the induction equa- 
tion in terms of the vector potential A, where B — V x A. 

We impose regularity conditions on the axis (vo — 0) and 
outflow boundary conditions on vo = n7 max and z — ±z max . 

2.2. A cool disc in a hot corona: The initial state 

A simple way to implement a cool, dense disc embedded in 
a hot, rarefied corona without modelling the detailed physics 
of coronal heating is to prescribe the specific entropy, s, such 
that s is smaller within the disc and larger in the corona. For a 
perfect gas this implies p = e s / ClJ q 1 (in a dimensionless form), 
with the polytrope parameter e s l Cv being a function of position. 
In the model considered here, we have an intermediate value for 
the specific entropy within the star. We prescribe the polytrope 
parameter to be unity in the corona and less than unity in the 
disc and in the star, so we put 

e s/c * = 1 - (l-/3disc)Cdisc ~ (1-/W)&tar, (5) 

where £disc and £ sta r are (time-independent) profiles specifying 
the shapes of the disc and the star. The free parameters, < 
/3disc, /?star < 1, control the entropy contrasts between disc and 
corona and between star and corona, respectively. 

In the absence of a magnetosphere, our initial state is a hy- 
drostatic equilibrium with no poloidal velocity, assuming an 
initially non-rotating hot corona that is supported by the pres- 
sure gradient. Since we model a disc that is cool, the disc is 
mainly centrifugally supported, and as a result it is rotating at 
slightly sub-Keplerian speed. 

The temperature ratio between disc and corona is roughly 
/3disc- Assuming pressure equilibrium between disc and corona, 
and p oc pT for a perfect gas, the corresponding density ratio is 
then (3~7r Thus, the entropy contrast between disc and corona 
chosen here {fidi SC = 0.005), leads to density and inverse tem- 
perature ratios of 200 : 1 between disc and corona. 

A rough estimate for the initial toroidal velocity, u„o, in 
the midplane of the disc follows from the hydrostatic equilib- 
rium as u v0 w y/1 — /SdiscMK- For /Jdisc = 0.005, the toroidal 
velocity is within 0.25% of the Keplerian speed. 

The initial hydrostatic solution is an unstable equilibrium 
because of the vertical shear between the disc, star and corona 
(Urpin & Brandenburg 1998). In addition, angular momentum 
transfer by viscous and magnetic stresses - the latter from the 
disc dynamo - drives the solution immediately away from the 
initial state. 

2.3. Application to a protostellar star-disc system 

Since we use dimensionless variables, our model can be 
rescaled and can therefore be applied to a range of different 
astrophysical objects. Here, we consider values for our nor- 
malization parameters that are typical of a protostellar star- 
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disc system. We scale the sound speed with a typical coro- 
nal sound speed of c s q = 10 2 kms -1 , which corresponds to 
a temperature of T w 4 x 10 5 K, and the disc surface den- 
sity with S = lgcm~ 2 . Further, we assume M* = 1M 
(where M* and M Q are the stellar and solar mass, respec- 
tively), a mean specific weight of /i = 0.6, and 7 = 5/3. 
This fixes the units of all quantities. The resulting velocity 
unit is [u] = c s o = 10 2 kms _1 , the length unit is [r] = 
GM*/[u} 2 m 0.1 AU, the time unit is [t] = [r]/[u] w 1.5 d, 
the unit for the kinematic viscosity and magnetic diffusivity is 
[v] = [rj\ = [u][r] ~ 1.5 x 10 19 cm 2 s _1 , the unit for specific 
entropy is [s] = c p = 7/(7-1)%/^ 3.5x 10 s cm 2 s~ 2 K _1 , 
the unit for specific enthalpy is [h] — [u] 2 — 10 4 km 2 s~ 2 , the 
temperature unit is [T] = [h]/[s] w 3 x 10 5 K, the density 
unit is [g] = So/M ~ 7.5 x 10~ 13 gcm~ 3 , the pressure unit 
is [p] = (7 — l)/7[p][/i] ~ 30gcm _1 s -2 , the unit for the 
mass accretion rate is [M] — So[w][r] ~ 2 x 1O~ 7 M yr -1 , 
the magnetic field unit is [B] = [it](47r[£>]) 1/2 w 30 G, and 
the unit for the magnetic vector potential is [A] = [B] [r] « 
4xl0 13 Gcm. 

Computations have been carried out in the domain (ro, z) € 
[0, 2] x [—1, 1], with the mesh sizes 6m = Sz = 0.01. In our 
units, this corresponds to the domain extending to ±0.1 AU in 
z and 0.2 AU in w. 

We choose a set of values for our model parameters such 
that the resulting dimensions of our system and the resulting 
initial profiles of the physical quantities are close to those for a 
standard accretion disc around a protostellar object. The stellar 
radius is = 0.15, corresponding to 3 solar radii, the disc 
inner radius is w- m — 0.6, i.e. 4 stellar radii, the disc outer 
radius is w ou t = 1.9, close to the outer domain boundary, and 
the disc semi-thickness is zq = 0.15, i.e. equal to the stellar 
radius. 

Further, our choice for the entropy contrast between disc 
and corona, /3disc = 0.005, leads to an initial disc temperature 
ranging between 9000 K in the inner part and 900 K in the outer 
part. Real protostellar discs have typical temperatures of about 
a few thousand Kelvin (e.g., Papaloizou & Terquem 1999). As 
turns out from our simulations, the disc temperature increases 
by less than a factor of 2 with time, except for the inner disc 
edge where the increase is much higher. The low disc temper- 
ature corresponds to a relatively high disc density. In the inner 
part of the disc, it is about 3 x 10~ 10 gem -3 initially, but in- 
creases to about 10~ 9 gcm~ 3 at later stages. In the outer part 
it is a few times 10~ n g cm~ 3 . 

The Shakura-Sunyaev coefficient of the turbulent disc vis- 
cosity is ass = 0.004 whereas the magnetic diffusivity is 
2 x 10~ 5 in the corona and in the star, and enhanced to 6 x 10 -5 
in the disc. The latter corresponds to ccgg = r)di S c/ (c s ,disc-Zo), 
ranging between 0.004 and 0.008 in the disc midplane. With 
?7disc = 6 x 10~ 5 (corresponding to about 9 x 10 14 cm 2 s _1 ) 
and the disc semi-thickness zq = 0.15, the diffusion time is 
tditt = z%/r/ disc = 375. 

The stellar surface angular velocity is about unity (in 
dimensionless units), corresponding to a rotation period of 
around 10 days. This means that the corotation radius is m co w 
1 in nondimensional units. The inner disc radius (which is fixed 
in our model) is rotating with roughly Keplerian speed, result- 



ing in a rotation period of the inner disc edge of around 4.5 days 
(3 in dimensionless units). 

2.4. Disc dynamo and stellar magnetosphere 

To date, almost all models of the formation and collimation of 
winds and jets from protostellar accretion discs rely on an ex- 
ternally imposed poloidal magnetic field and ignore any field 
produced in the disc. This is not the case in the model devel- 
oped in Paper I, which also forms the basis of the model used 
in the present paper. In Paper I we study outflows in connection 
with magnetic fields that are solely generated and maintained 
by a disc dynamo, resolving the disc and the star at the same 
time as well as assuming a non-ideal corona. There, the ini- 
tial seed magnetic field is large scale, poloidal, of mixed parity, 
weak and confined to the disc. 

However, observations of T Tauri stars suggest that mag- 
netic star-disc coupling might result in a spin-down of the star 
due to magnetic braking by a stellar magnetosphere penetrating 
the disc. Magnetic fields as strong as 1 kG and larger have been 
detected on T Tauri stars (e.g., Guenther et al. 1999), and there 
is evidence for hot and cool spots on the stellar surface. 

Our present model is similar to that of Paper I except that, in 
addition to the disc dynamo, we also model a magnetosphere of 
the protostar. For the magnetosphere, we assume that in the ini- 
tial state a stellar dipolar magnetic field threads the surrounding 
disc, leading to a contribution of A v given by 



A v {zu, z) = ... + Ajtar— Y~ (1 - fstar) , (6) 

where r* is the stellar radius, r = (w 2 + z 2 ) 1 / 2 is the spheri- 
cal radius, and A star is a parameter controlling the strength of 
the stellar magnetosphere. Since m = r sin 9 (with 9 the co- 
latitude), the magnetic moment is /i mag = A stav r 2 (cf. Model I 
in Miller & Stone 1997). The dots in Eq. (|6} denote the con- 
tribution from the seed magnetic field in the disc. Since the 
magnetic field describing the magnetosphere is force-free, the 
hydrostatic equilibrium is not affected but equal to the magne- 
tostatic equilibrium. £ st ar is a smoothed profile for the star and 
equal to unity only at the origin (centre of the star), so that the 
magnetosphere extends to parts of the star as well. In order to 
anchor the magnetosphere within the star and on the stellar sur- 
face, q, u, and A are set to their initial values at each timestep 
in the region including the star and extending to about 1.5?v 
Therefore, in the models with magnetosphere, there is no mass 
sink in the star of the type discussed in Paper I. The star is re- 
solved in that it stays within the mesh, allowing the study of the 
effects of a spinning dipole (cf. Sect. 12. 1> . Note that outside the 
anchoring region, the magnetosphere is not imposed and can 
evolve dynamically with time. In the induction equation, we do 
not include any term containing an external magnetic field. As 
a reference model we choose A stal — 5, which corresponds 
to a stellar surface magnetic field strength of about 1 kG. The 
magnetic moment is then about 10 37 G cm 3 . 
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Fig. 2. Poloidal velocity vectors and poloidal magnetic field 
lines (white) superimposed on a colour/grey scale representa- 
tion of log 10 h for our reference model (Model Ml) with a stel- 
lar surface magnetic field strength of about 1 kG and with disc 
dynamo at the time t = 1267. Specific enthalpy h is directly 
proportional to temperature T, and log 10 h = (—2, — 1, 0, 1) 
corresponds to T w (3x 10 3 , 3x 10 4 , 3x 10 5 , 3x 10 6 ) K. The 
black dashed line shows the surface where the poloidal veloc- 
ity equals the Alfven speed from the poloidal magnetic field 
( Alfven surface). The disc boundary is shown as a thin black 
line. In red are marked the stellar surface as well as the surface 
up to which the magnetosphere is anchored. Here and in the 
following figures, the field lines tend to become vertical near 
the boundaries, but this is due to artifacts from the boundary 
conditions. 

3. Results 

We consider models with magnetospheres of different 
strengths. We declare as our reference model one with a stel- 
lar surface field strength of 1 kG (Sect. 13. 1> . and compare it 
with a model without disc dynamo (Sect. 13. 21 . and with mod- 
els whose stellar surface field strength is varied between zero 
and 5kG (Sects [33M3.7I . Finally, the magnetic and accretion 
torques on the star are studied for different magnetospheric 
models (Sect.l3~8l. 

3.1. Reference model: \B Bmf \ r* lkG (Model M1) 

Our simulations show that, similar to Paper I, also the present 
model with a stellar magnetosphere (in addition to the disc 
dynamo) develops a structured outflow, composed of a stel- 
lar wind and a disc wind; see Fig. |2] The disc wind is faster, 
cooler and less dense, with the highest velocities in the wind 
originating from the inner edge of the disc, whereas the stel- 
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Fig. 3. Reference model (Model Ml) at the time t = 1267 
when the accretion rate is maximum. Colour/grey scale repre- 
sentation of the density (bright colours or light shades indicate 
high values; dark colours or dark shades indicate low values) 
with poloidal magnetic field lines superimposed (white) and 
the azimuthally integrated mass flux density, represented as the 
vector 2ir-cag(u-cj, u z ), shown with arrows (except in the disc 
where the density is high and the mass flux vectors would be 
too long). 



lar wind is slower, hotter and denser. The inner disc wind is 
magneto-centrifugally accelerated, whereas the outer disc wind 
and the stellar wind are mostly pressure-driven. These driving 
and acceleration mechanisms will be explained in Sect. 13. 51 us- 
ing Model S in which the outflow structure is most pronounced. 

The disc wind velocity reaches about 240 kms^ 1 , whereas 
the stellar wind velocity goes only up to about 10 kms -1 , The 
disc wind mass loss rate is time-dependent with an average 
value of around 2 x 10 -7 Mq yr" 1 . 

Accretion of matter onto the central star is highly episodic 
with a maximum rate of about 4 x 10 _9 M Q yr _1 . Figure [3] 
shows a snapshot at a time when the accretion rate is max- 
imum. We estimate the mass accretion rate in the following 
way. We put a cylinder around the axis extending to vj = 0.5 
and z — ±0.3. (Note that the disc inner edge is at m = 0.6, 
and the magnetosphere is anchored up to ro w 0.275 due to 
a smoothed profile for the anchoring region.) Then we calcu- 



6 



B. von Rekowski and A. Brandenburg: Outflows and accretion in a star-disc system 




m 0,0 "I 



. 1 ; n 

0.0 


0,5 


1.0 


1.5 2,0 






a 







2.5 2.0 1.5 1.0 5 0.0 0.5 1.0 



Fig. 4. Model Ml-0. Same model as the reference model in 
Fig. 13 but with no disc dynamo, i.e. a = in Eq. 0, and 
averaged over times t = 1100 . . A = 1127. Note that, unlike 
all the models with disc dynamo, the magnetic field (field lines 
shown in white) is almost entirely swept out of the disc. 



late the azimuthally integrated mass flux density, 2-kzdqu, nor- 
mal to the boundaries of the cylinder and integrate it over the 
boundaries. We take into account that there is mass loss due 
to the stellar wind. It turns out that in this model, matter en- 
ters the cylinder mainly through the vertical boundary and only 
then partly flows along magnetospheric field lines, because the 
magnetosphere does not extend far enough. 

3.2. A model with no disc dynamo (Model M1-0) 

Almost all previous work on the magnetic star-disc coupling 
assumes that the magnetic field in the disc results entirely from 
the central star. Some models include an externally imposed 
magnetic field (e.g., Kiiker et al. 2003). In order to study the 
effect of the disc dynamo on the overall field structure and on 
the resulting outflow, we have calculated a model where the 
disc dynamo is turned off, i.e. we put a = in Eq. l|3}. 

As can be seen in Fig. @] the magnetic field in the disc 
is now almost entirely swept away. The initial disc field due 
to the penetrating magnetosphere is expelled and only a very 
weak disc field remains without a dynamo. As a consequence, 
the structure and driving/acceleration mechanism of the disc 
wind change significantly. Transport of specific angular mo- 
mentum is weak and almost entirely in the w direction, along 
very weak field lines threading the disc and running in the low 
corona almost parallel to the disc surface. The resulting disc 
wind along these field lines is therefore also mainly in the hor- 
izontal direction (Fig. 0}. Since the field threading the disc is 



very weak, the disc wind is always super-Alfvenic, reaching 
a speed of about 160 kms -1 , so that no magneto-centrifugal 
acceleration can take place. Characteristic for the models with 
disc dynamo is a cooler, less dense region ("conical shell") that 
originates from the inner disc edge, and that has higher spe- 
cific angular momentum than elsewhere and contains the inner 
magneto-centrifugally accelerated disc wind. This conical shell 
disappears as the magnetic field, that was threading the disc, is 
swept away and no disc dynamo is generating and maintaining 
a strong enough field. As a result, the pressure force becomes 
then the only driving mechanism of the whole disc wind. 

As the disc field is expelled, the size of the closed stellar 
magnetosphere shrinks; the extended open stellar field lines 
now find place to run in the corona parallel above the very 
weak field lines threading the disc; see Fig. |4] [Note that in 
our magnetospheric models, the vertical alignment of the field 
lines near the boundaries is due to boundary conditions. In or- 
der to stabilize our code at the boundaries (especially in situa- 
tions where the wind speed is low), we force matter to be slowly 
advected from the computational domain normal to the bound- 
aries, rather than in spherical radial direction as in Paper I.] The 
stellar wind, however, does not seem to be much affected, with 
a maximum speed of about 10 kms -1 . 

The average disc wind mass loss rate does not change much 
in comparison with the reference model, and remains around 
2 x 1O~ 7 M yr _1 . Accretion is still highly episodic; the max- 
imum mass accretion rate is roughly 5 times larger compared 
to the same model but with disc dynamo, so it is now about 
2 X 10- 8 M Q yr- 1 . 



3.3. Weak magnetosphere: 
(Model W) 



\B smi \ «200G 



A stellar surface field strength of 200 G is not large enough 
to maintain a closed magnetosphere; see Fig. [5] The field 
lines open up entirely to form open stellar and disc fields. 
This mainly affects the stellar wind; its maximum speed is 
reduced to about 2.5kms~ 1 , whereas the disc wind velocity 
reaches about 200 kms -1 . However, the disc dynamo again 
produces a structure in the disc wind with different driving 
mechanisms. Transport of specific angular momentum is also 
enhanced along field lines threading the midplane at a radius 
smaller than the prescribed inner disc radius, which is due to 
an increased angular velocity. In the disc, the field lines show 
a characteristic wiggly pattern similar to the so-called chan- 
nel flow solution seen in two-dimensional simulations of the 
Balbus-Hawley instability (Hawley & Balbus 1991). This in- 
stability is indeed to be expected for the weak field strengths 
in the present model; for stronger fields this instability is sup- 
pressed. Note, however, that this channel flow type solution in 
the disc can only be seen in the time-averaged picture shown in 
Fig.0 

In this connection it might be worthwhile reiterating that 
the purpose of introducing a turbulent viscosity is of course 
to imitate the effects of the three-dimensional Balbus-Hawley 
instability. The channel flow solution is clearly an artifact of 
two dimensions, and we regard the almost complete absence of 
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Fig. 5. Model W (weak magnetosphere and disc dynamo); the 
magnetic field strength at the stellar surface is about 200 G. 
Averaged over times t — 825 . . . t — 850. Shown is the same 
as in Fig- EI 
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Fig. 6. Model M2 (medium magnetosphere and disc dynamo); 
the magnetic field strength at the stellar surface is about 2 kG. 
t = 562. Shown is the same as in Fig. EJ 

the channel flow behaviour in our solutions as a confirmation 
that there is no "double-instability" of the mean flow in the 
sense discussed above [see also Tuominen et al. (1994) and 




0.0 0.2 0.4 0.6 0.8 1.0 

CD 



Fig. 7. Model M2. t = 562. Shown is the same as in Fig. [3] 
also at a time when the accretion rate is maximum. 



Yousef et al. (2003) for discussions regarding the problems of 
stability studies of mean-field models]. 

3.4. Medium magnetosphere: \B smi \ rj 2kG 
(Model M2) 

A stellar magnetosphere with an increased stellar surface field 
strength of about 2 kG interacts with the disc dynamo in such a 
way that the resulting magnetic field, temperature and density 
distributions are clearly influenced. The structure in the disc 
wind becomes more visible, with a more pronounced (cooler 
and less dense) conical shell originating from the inner disc 
edge containing a (slightly) faster magneto-centrifugally ac- 
celerated inner disc wind with a maximum speed of about 
250 kms -1 ; compare Figs[2]and[6] The average disc wind mass 
loss rate stays around 2 x 10~ 7 A/q yr _1 , and the stellar wind 
velocity does not exceed 10 km s _1 . 

When the stellar surface field is as strong as 2kG, the 
highly episodic accretion is correlated with magnetospheric os- 
cillations. During time intervals when the magnetosphere is ex- 
panded so that the outer closed field lines are sufficiently close 
to the disc, the accretion flow is along magnetospheric field 
lines; see Fig. (see also Sect. I3.5l l. The maximum accretion 
rate is about 2 x 10~ 8 Af© yr _1 . 
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Fig. 8. Model S (strong magnetosphere and disc dynamo); the 
magnetic field strength at the stellar surface is about 5kG. 
Mass accretion rate estimated as described in Sect. 13.11 (up- 
per panel) and vertical magnetic field component at (w, z) = 
(0.6, 0.4) (lower panel) as functions of time. The latter is an 
indicator for the oscillations of the magnetosphere. Note that 
the accretion rate and the magnetospheric oscillations are cor- 
related such that the accretion rate is high when B z at the posi- 
tion given above is low. This corresponds to the state when the 
star is connected to the disc by its magnetosphere. Time is here 
in non-dimensional units, but with a time unit of 1 .5 days, the 
oscillation period is seen to be around 15 to 30 days. 
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3.5. Strong magnetosphere: 
(Model S) 
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The time sequence shown in the lower panel of Fig. [8] shows 
that the magnetosphere is oscillating with a period of around 
15 to 30 days. As we shall show below, the oscillations mean 
that the magnetosphere is expanding and contracting and in this 
way periodically connecting and disconnecting open stellar and 
disc fields, therefore changing the magnetic star-disc coupling. 

It turns out that the highly episodic accretion flow and ac- 
cretion rate are correlated with the configuration of the mag- 
netosphere (compare upper and lower panels of Fig. [8}. Low 
values of B z at the position (ru, z) — (0.6, 0.4), correspond 
to a configuration where closed magnetospheric lines penetrate 
the inner disc edge, thus connecting the star to the disc (as in 
Fig- ED- In this case, disc matter is loaded onto magnetospheric 
field lines and flows along them to accrete onto the star. As a 
consequence, the accretion rate is highest during these phases, 
typically up to between 10~ 8 Af© yr _1 and 2.5x 10~ 8 Af© yr _1 
(see Fig. [8] upper panel). 

High values of B z at the position (zu, z) = (0.6, 0.4), cor- 
respond to a configuration when the outer field lines of the mag- 
netosphere have opened up into disconnected open stellar and 
disc field lines, thus disconnecting the star from the disc (as in 
Fig.llOi. In this case, matter is lost directly into the outflow and 
there is no net accretion of disc matter (see Fig. [8] upper panel). 

However, even when the accretion rate is high, most of the 
mass that leaves the disc goes into the wind. The mass loss 
rates into the disc wind are about one order of magnitude higher 



Fig. 9. Model S. Shown is the same as in Fig. [3] but at the time 
t = 98 when the star is magnetically connected to the disc. 



compared to the peak accretion rates (see Fig.[^ upper panel), 
and not very different from the rates in our other magneto- 
spheric models presented in previous subsections, ranging be- 
tween 10" 7 M© yr^ 1 and 2.5 x 10~ 7 Af© yr _1 . The disc wind 
mass loss rates are fluctuating with time on rather short time 
scales, and no correlation with the oscillating magnetosphere 
can be detected (compare upper and lower panels of Fig. II 1>. 

The time-dependent behaviour of the accretion flow and the 
disc wind is interesting in view of observations. Recent high- 
resolution short- and long-term monitoring of classical T Tauri 
stars (CTTS) with the ESO Very Large Telescope reveals that 
CTTS are highly variable on time scales of minutes to several 
years. This variability in the emission spectra is associated with 
both accretion and outflow processes (Stempels & Piskunov 
2002, 2003). 

As in all our models with disc dynamo, the dynamo pro- 
duces a structured disc wind. However, the larger the stellar 
surface magnetospheric field strength is, the more pronounced 
is the structure. The outflow consists of (i) a slower, hotter and 
denser, mostly pressure-driven stellar wind, (ii) a faster, cooler 
and less dense outer disc wind that is also mostly pressure- 
driven, and (iii) a fast magneto-centrifugally accelerated inner 
disc wind within a clearly visible conical shell originating from 
the inner disc edge, that has lower temperatures and densities 
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Fig. 10. Model S. Shown is the same as in Fig. [9] but at the 
time £ = 102 when the star is magnetically disconnected from 
the disc. 

than elsewhere (see Fig.ll2>. This outflow structure was already 
discussed in Paper I in some detail. 

Figure^]shows a snapshot at a time in a transition period, 
when the magnetic star-disc connectivity changes. In these pe- 
riods when disconnected stellar and disc field lines are about to 
reconnect, in addition to the above described structured outflow 
typical for our models with disc dynamo, there is a pressure- 
driven, hot and dense, but relatively fast outflow between the 
stellar and disc winds. 

As one can see in Figs [D] and ^] specific angular mo- 
mentum is carried outwards from the disc into the corona 
mainly along magnetic field lines threading the innermost part 
of the disc. Here the dynamo-generated field is strongest and 
is responsible for confining the conical shell in the corona. 
The angle between the rotation axis and these field lines is 
30° and larger at the disc surface, which is favourable for 
magneto-centrifugal acceleration (Blandford & Payne 1982; 
see Campbell 1999, 2000, 2001 for a more detailed treatment). 
Furthermore, in the conical shell the flow is highly supersonic 
but the Alfven radius is about two times larger than the radius 
where the field lines have their footpoints at the disc surface. 
This is sufficient for the magnetic field lines to act as a lever 
arm along which the inner disc wind is accelerated. On the 
other hand, in the outer disc wind the Alfven surface is very 
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Fig. 11. Model S. Disc wind mass loss rate (upper panel) and 
vertical magnetic field component at (vj,z) = (0.6, 0.4) (lower 
panel) as functions of time. The latter is an indicator for the 
oscillations of the magnetosphere. No correlation between the 
disc wind mass loss rate and the magnetospheric oscillations 
can be seen. 

close to the disc surface. As a consequence, acceleration due to 
the pressure gradient becomes more important there, although 
the criterion of Blandford & Payne (1982) is still fulfilled. The 
stellar wind is mostly pressure-driven, carrying almost no spe- 
cific angular momentum and remaining subsonic. The outflow 
between the stellar and disc winds has small azimuthal veloc- 
ity and small specific angular momentum, and becomes quickly 
super- Alfvenic. 

A confirmation of the presence of coexisting pressure- 
driving and magneto-centrifugal acceleration of the winds as 
described above, can be obtained by looking at the ratio be- 
tween the poloidal magneto-centrifugal and pressure forces, 



,(p)i 



I j. j |/ |x ■ poll' wnere the subscript 'pol' denotes the poloidal 
components. The forces appear in the equation of motion as 



g(n 2 zv-V$) +JxB, 



(7) 



and F^ p ' = — Vp. Their ratio is shown in Fig. [21 Again, the 
conical shell can be clearly seen. In the conical shell, the ratio 
assumes large values which confirms that magneto-centrifugal 
launching and acceleration of the disc wind is dominant there. 
On the other hand, at the outer parts of the disc surface the pres- 
sure force is stronger, leading to pressure-driving. At the stel- 
lar surface, both forces have comparable strengths; however, in 
large parts of the stellar wind the pressure force becomes dom- 
inant. The outflow between the stellar and disc winds is also 
pressure-dominated. 

Figures^]andE]show the distribution of specific angular 
momentum for the two cases, where the star is connected to and 
disconnected from the disc, respectively. Comparing these two 
figures, one can see that in periods of no net accretion, when 
the star is disconnected from the disc, specific angular momen- 
tum transport by the inner disc wind is visibly enhanced, car- 
ried by enhanced stellar and inner disc winds. Their velocities 
are larger in these periods: terminal outflow speeds are then 
about 20kms~ 1 in the stellar wind and about 400 kms -1 in 
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Fig. 12. Model S. t — 103 when the star is about to be magnet- 
ically reconnected to the disc. Shown is the same as in Fig. 13 




Fig. 13. Colour/grey scale representation of the specific angu- 
lar momentum, I = wu v — w 2 VL, normalized by the max- 
imum specific angular momentum in the disc, Zdi sc ,max, with 
poloidal magnetic field lines superimposed (white) for Model S 
at the time t = 98 when the star is magnetically connected to 
the disc. The blue solid line shows the Alfven surface, and the 
orange dashed line the sonic surface. 



Fig. 14. Same as in Fig. \H\ (Model S), but at the time t = 102 
when the star is magnetically disconnected from the disc. Note 
that now a lot of specific angular momentum leaves the disc 
from its inner edge. 



the conical shell of the disc wind, where velocities are highest. 
Accretion flow velocities in the disc are up to about 15 kms -1 . 
In our computational domain, both the inner stellar wind and 
the funnel flow remain subsonic at all times. The slow stellar 
wind is likely to be due to our boundary conditions that we 
impose in our magnetospheric models by fixing the poloidal 
velocity to be zero in the star. 

A peculiar feature of all our models with a disc as cool as 
in the models of this paper, is the presence of at least one radial 
polarity reversal of the dynamo-generated magnetic field in the 
disc (cf. Figs 23-25 of Paper I). In the model with strong mag- 
netosphere discussed in this section, the first reversal occurs 
roughly at the corotation radius (m co w 1), where the angular 
velocity is about 7.5 x 10 _6 s , corresponding to a rotation 
period of around 10 days. 

The results in Sects 13.41 and 13.51 concerning the oscillat- 
ing magnetosphere and associated reconnection processes and 
episodic accretion, look indeed like being in the regime de- 
scribed by Goodson & Winglee (1999). However, our stellar 
wind is relatively slow and mostly pressure-driven whereas 
their stellar wind (jet) is fast and magneto-centrifugally driven. 
On the other hand, the main new feature induced by our 
disc dynamo is the structure in our disc wind. Whereas the 
outer disc wind is slower and mostly pressure-driven, the in- 
ner disc wind is fast, and launched and accelerated magneto- 
centrifugally; it might collimate at larger heights to form the 
observed protostellar jet. The magnetospheric extension is lim- 
ited by the dynamo-generated disc field that is strongest in the 
inner disc parts and advected into the disc corona. 
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Fig. 15. Colour/grey scale representation of the ratio be- 
tween the poloidal magneto-centrifugal and pressure forces, 

l^poi^l/l-Fpoll as defined in Eq. Q and below, with larger 
values corresponding to lighter shades. Superimposed are the 
poloidal magnetic field lines (white). The blue solid line shows 
the Alfven surface. Model S at the time t = 103 when the star 
is about to be magnetically reconnected to the disc. 
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Fig. 16. Poloidal magnetic field vectors of the model shown 
in Fig. El Note that the magnetic field reversal in the corona 
indicates the presence of a current sheet between the magneto- 
sphere and the field threading the disc at about 45°. This means 
that no X-point forms. The magnetic field vectors are not shown 
in the anchoring region. The length of the vectors is weighed 
with m. The black dashed line shows the Alfven surface. 



3.6. Dependence on stellar field strength 



Table 1. Summary of parameter values related to the physics 
of both the magnetosphere and the disc dynamo for a sequence 
of models where the strength of the stellar field is varied. The 
hyphen in the row for M 1-0 indicates that there is no turbulence 
in that run. Note that |_B S urf | = Aitar/r* x [B]. 
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In Table ^ we summarize the parameters that have been 
changed in the different models considered in this paper and 
that are related to the physics of both the magnetosphere and 
the disc dynamo. The strength of the stellar field is given both 
in nondimensional and in dimensional units, together with the 
ao coefficient quantifying the strength of the a effect (ao < 
means that a is negative in the upper disc plane). 



As the strength of the stellar field is reduced, the size of the 
stellar magnetosphere shrinks. At the same time, the field in 
the disc also decreases somewhat. Although the conical struc- 
ture becomes more pronounced as the strength of the stellar 
field is increased, the overall magnetic field and outflow struc- 
tures remain the same for medium and strong stellar surface 
magnetospheric field strengths. In particular, the field struc- 
ture (inside the corotation radius) is not of X-point topology, in 
contrast to what is described by Shu et al. (1994). Instead, the 
dynamo-generated magnetic field in the disc always arranges 
itself such that, in the corona, the field threading the disc is 
anti-aligned with the stellar dipole. This necessarily leads to 
magnetospheric current sheets; see Fig. 1161 

3.7. A model with no magnetosphere (Model N) 

In this section we present a model without stellar magneto- 
sphere where the only source of magnetic fields is the dynamo 
operating in the disc. The field threading the star results en- 
tirely from advection of the dynamo-generated disc field and is 
therefore maintained by the disc dynamo. In this model without 
magnetosphere, there are no restrictions to the magnetic field in 
the star. Figure ^] shows a time-averaged picture of the same 
model as in Fig. 25 of Paper I, but at later times, when the disc 
dynamo is saturated. Note that in this model the inner and outer 
disc radii are smaller than in the models with magnetosphere. 
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Fig. 17. Similar to Fig. [2] but for the same model as in Fig. 25 
of Paper I, where now we have averaged over later times 
t = 584 . . . t = 590, when the disc dynamo is saturated. This 
corresponds to Model N (no magnetosphere and disc dynamo) 
in this paper. (In this figure, the black dashed line shows the 
surface where the poloidal velocity equals (c s 2 + v\ poi) 1 ^ 2 , 
with c s the sound speed and UA,poJ the Alfven speed from the 
poloidal magnetic field.) Note the similarity of the overall out- 
flow and coronal magnetic field structures to Model S; Fig. 1121 



The overall outflow and coronal magnetic field structures 
are very similar to Model S fFig. H2l> . as are the driving mecha- 
nisms, cf. Fig.^Dwith Figs^Jand^] However, as one can see 
in Figsll7landll8l here the stellar wind is much faster and be- 
comes supersonic; in contrast to Model S, in this model without 
magnetosphere the velocity can freely evolve in the star, and 
the only prescription for the star is to act as a self-regulatory 
mass sink (cf. Sect. 12. li . The terminal stellar wind velocity is 
about 150 kms -1 , and the terminal inner disc wind velocity in 
the conical shell is also higher (about 500 kms -1 ). 

Accretion is not episodic as in the magnetospheric 
models but the accretion rate is always positive (around 
1O _6 M0 yr _1 ), and the accretion flow is in the cylindrical 
radial direction. The disc wind mass loss rate is around 4 x 
1O _7 M0 yr -1 . Both rates are therefore higher than the rates 
in the models with magnetosphere. The disc wind mass loss 
rate is only about twice as large, but the accretion rate is at 
least 50 times larger, compared to when the episodic accretion 
in the magnetospheric models reaches its maximum rate. This 
changes the ratio between mass accretion rate and disc wind 
mass loss rate not only quantitatively, but also qualitatively. 
This ratio is now about 7 : 3, but it was less than its reciprocal 
value (about 1 : 9 at most) in the models with magnetosphere. 



Fig. 18. Shown is the same as in Figs [O] and ^] but for the 
same model (Model N) and time as in Fig. 1171 



Nevertheless, the fact that the stellar magnetic field is still 
able to shield a relatively large part of the accretion flow and 
deflect it into the wind could be the result of the assumption of 
fully axisymmetric flows. Actual T Tauri stars possess highly 
nonaxisymmetric, non-dipolar fields (Johns-Krull et al. 1999, 
Johns-Krull & Valenti 2001). However, the ratio of accretion 
rate to wind mass loss rate is difficult to determine observa- 
tionally with great confidence. Although Pelletier & Pudritz 
(1992) estimate that only about 10% of the matter joins the 
wind, our ratios lie still in the range of observationally derived 
ratios, even those for our models with magnetosphere. 

3.8. Magnetic and accretion torques 

In an axisymmetric system, the equation of conservation of an- 
gular momentum reads 

d 

-q^KQI) = -V ' (*acc + *mag + *visc), (8) 

where I — vou v — ru 2 Q, is the specific angular momentum, 

t acc = £>G7WpollV = QU vo \W 2 n (9) 

is the material stress (it po i is the poloidal velocity field), 

, -Bpol-By ,..„, 

tmag = (10) 

47T 

is the magnetic stress (B po \ is the poloidal magnetic field), and 

*visc = -£^ t TO 2 VO (11) 
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Fig. 19. Dependence of the net accretion torque T acc net (upper 
panel) and the magnetic torque T mag (middle panel) on time at 
the spherical radius r « 0.33 (close to the star) for Model S. 
Lower panel: same indicator for the magnetospheric oscilla- 
tions as in Fig. [H] lower panel, and Fig.[^ lower panel. 

is the viscous stress (z/ t is the turbulent kinematic viscosity). 
At any given time, the total torque T acting on a central object 
(star) is given by the volume integral 



T 



jv- { t 



acc ~r ^mag ^visc 



■dV, 



(12) 



where the integral has to be taken over a volume enclosing the 
star. This is equivalent to the surface integral 



T = — f(t 



3 ).dS 



(13) 



= j{-QU vo ivj 2 n + vo Bpo ^ v + QV t VJ 2 Vtt) ■ dS, (14) 

where dS is the outward directed surface element. In other 
words, the total torque acting on a star is given by the angular 
momentum flux across a surface enclosing the star (see Ghosh 
& Lamb 1979b). If the flux is negative, i.e. towards the star, 
then the torque on the star is positive. 

We calculate the angular momentum flux towards the star - 
carried by the matter (T acc ), by the magnetic field (T mag ), and 
due to the viscous stress (T v i sc ) - across spheres around the star, 
so that we can write in spherical polar coordinates (r, 9, (p) 
(with 9 the co-latitude, w = r sin 9 and e r the radial unit 
vector): 



T acc (r) = -2vrr 2 f (t acc • e r ) sin9 d9 
Jo 



-27rr / gu r m fi sin 9 d9 



for the accretion torque, 



TmagM = -27rr 2 / (t mag • e r ) sin 9 d9 
Jo 



2irr~ I vo 

lo 4^ 



BrBv sin9d9 
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Fig. 20. Magnetic torque T mag as function of spherical radius 
r for the magnetospheric models presented in this paper. T mag 
is shown for radii between the approximate anchoring radius 
ro = 0.275 and the inner disc radius zu = 0.6. Purple dash- 
triple-dotted: Model W, averaged over times t = 825 . . . t — 
850; orange dashed: Model Ml, t = 1267; green dash-dotted: 
Model Ml-0, averaged over times t = 1100 . . . t — 1127; blue 
long-dashed: Model M2, t = 562; red solid: Model S, t = 
103. The black dotted line marks the magnetic torqueless state 
(T mag = 0). 



for the magnetic torque, and 

T visc (r) = ~2nr 2 / (i visc • e r ) sin 9 d9 
Jo 

= 2irr I QVtvo 2 -— sin9d9 
Jo or 



(19) 
(20) 
1 corre- 



for the viscous torque. In our units, T[ acc mag v ; sc i 
sponds to about 2.5 x 10 38 erg. 

At each position between the star and the disc where mag- 
netic field lines are connected to the star, B r B v > means 
that those field lines are leading the star in a rotational sense, 
therefore spinning the star up. Therefore, T mag > translates 
into a magnetic spin-up of the star. Conversely, B r B v < 
means that those field lines are lagging behind the star and are 
therefore spinning the star down, so T mag < translates into a 
magnetic spin-down of the star. 

However, — u r voVl > means only that angular momen- 
tum is added to the star. Whether or not this also means a stellar 
spin-up (or spin-down if — u r wQ, < 0), depends on the magni- 
tude of the angular velocity of the matter, f2(r, 9), compared to 
the magnitude of the effective rotation rate of the stellar surface, 
f2* (see AppendixlAlfor a precise definition). We therefore de- 
fine the net accretion torque, 



t(r) 



PIT 

-2vrr 2 / gu r zu 2 {fl-n*) sin OdQ. (21) 
Jo 



If the accreting matter (u r < 0) is rotating faster than the star, 
the star will be spun up and T ac c,net > 0. Likewise, if the ac- 
creting matter is rotating slower than the star, the star will be 
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Fig. 22. Net accretion torque (upper panel) and magnetic torque (lower panel) as functions of latitude at the spherical radius 
r « 0.33 (close to the star) for Model S for the two distinctive states: when the star is magnetically connected to the disc (t = 98, 
left) and when the star is magnetically disconnected from the disc (t = 102, right). Latitudes marked with "acc" are regions of 
strongest mass accretion (cf. Fig. [9} and latitudes marked with "out" are regions of strongest mass outflow (cf. Fig. 1101 ). 



magnetic torque 




0.30 0.35 0.40 0.45 0.50 0.55 0.60 
spherical radius 

Fig. 21. Magnetic torque T mag as function of spherical radius 
r for Model S at three different times. T mag is shown for radii 
between the approximate anchoring radius tq = 0.275 and the 
inner disc radius w = 0.6. Blue long-dashed: t — 98 when 
the star is magnetically connected to the disc; orange dashed: 
t = 102 when the star is magnetically disconnected from the 
disc; red solid: £ = 103 when the star is about to be magneti- 
cally reconnected to the disc (cf. Fig.l20>. The black dotted line 
marks the magnetic torqueless state (T mag = 0). 



spun down and T acc no t < 0. On the other hand, a stellar wind 
(u r > 0) that is rotating faster than the star, will spin down 
the star, and T acc . no t < 0, while a stellar wind that is rotating 



slower than the star, will spin up the star, so T) 



acc, net 



> 0. 



Figure [T9l shows that in Model S, close to the star where 
magnetic field lines are connected to the star, the net accretion 
torque is fluctuating around zero, whereas the magnetic torque 
is positive at all times. The magnetic torque is always much 
larger than the sum of the net accretion and viscous torques 



(T v i sc is not shown), suggesting a total spin-up of the star, the 
angular momentum flux towards the star being mainly carried 
by the magnetic field. A correlation is clearly visible between 
Tmag and the magnetospheric oscillations, and therefore also 
between T mag and the mass accretion rate. When the star is 
connected to the disc by its magnetosphere, i.e. when B z is 
lowest, both the accretion rate and the magnetic spin-up torque 
are maximum. 

The peak values of the magnetic spin-up torque are around 
2 x 10 39 erg, while those of the material spin-up torque are 
around 2.5 x 10 37 erg, and those of the material spin-down 
torque are around 5 x 10 37 erg. 

Also Romanova et al. (2002) find that most of the angular 
momentum flux to the star is carried by the magnetosphere. In 
the case of a protostar rotating with a period of around 9.4 days 
and a stellar surface field strength of about 1.1 kG, they find 
that matter carries only about 1% of the total flux. The mag- 
netic torque is positive so that it acts to spin up the star, and its 
amplitude is also correlated with the accretion rate. 

In all our models with sufficiently strong magnetosphere 
(|-B surf | > 1 kG; Models Ml, M1-0,M2 and S),T mag is mainly 
positive between the anchoring region and the inner disc edge 
(0.275 < zu < 0.6); see Fig.|20| This is roughly the region 
where magnetic lines are connected to the star, so that the 
star experiences the torque. This indicates that for sufficiently 
strong magnetosphere, angular momentum will be added to the 
star, carried by the magnetic field, resulting in a stellar spin-up 
by the magnetic field. As the stellar surface field strength of 
the magnetosphere increases, also the magnetic torque at radii 
between star and disc generally increases. In this region, the 
accretion and viscous torques (not shown) are negligible com- 
pared to the magnetic torque in all our magnetospheric models. 

Figure |^ shows that when star and disc are magnetically 
connected, the magnetic spin-up torque is higher at all spheri- 
cal radii between the approximate anchoring radius ro = 0.275 
and the inner disc radius w = 0.6. Basically what happens is 
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this: when the star is connected to the disc, there is accretion 
along magnetospheric field lines, and a large positive magnetic 
torque leads to a magnetic spin-up of the star. Then the star dis- 
connects from the disc, and the inner disc edge is magnetically 
spun up (see Fig. l2U : an enhanced inner disc wind carries away 
excess specific angular momentum (see also Fig. ll4t . 

The latitudinal dependencies of the torques for Model S, 
shown in Fig. [22] reveal that there is also a correlation be- 
tween the net accretion torque [i.e. the non-integrated quantity 
of Eq. J21H and the mass accretion rate. At latitudes of accre- 
tion onto the star, the net accretion torque is negative (i.e. a 
material spin-down), whereas at latitudes of stellar outflow, the 
net accretion torque is positive (i.e. a material spin-Mp). This is 
due to the small rotation rates between the star and the disc. 

4. Conclusions 

In agreement with earlier work by Hirose et al. (1997), 
Goodson & Winglee (1999) and Matt et al. (2002), our work 
confirms the possibility of episodic and recurrent magneto- 
spheric accretion also for models where the disc magnetic 
field is not imposed but dynamo-generated. The critical stel- 
lar surface field strength required for the episodic accretion to 
be correlated with the magnetic star-disc coupling, is around 
(or below) 2 kG, which is well in the range of observed field 
strengths for T Tauri stars (Johns-Krull et al. 1999). These re- 
current changes in the connectivity between the stellar magne- 
tosphere and the dynamo-generated disc field result in episodic 
mass transfer from the disc to the star. For the same stellar sur- 
face field strengths, the time-dependent wind velocities are also 
correlated with the magnetic star-disc coupling. The wind ve- 
locities as well as specific angular momentum transport from 
the disc inner edge into the corona are enhanced during peri- 
ods when the star is disconnected from the disc. Highly time- 
dependent accretion and outflows have also been detected ob- 
servationally in classical T Tauri stars (CTTS) by Stempels & 
Piskunov (2002, 2003). 

The observed stellar fields are, however, not dipolar, but 
show a strong nonaxisymmetric component, although the ge- 
ometry of the magnetosphere of protostars is yet unknown. 
[The assumption of a stellar dipolar field is motivated by ob- 
servations that suggest that the stellar field might be concen- 
trated at the poles in rapidly rotating CTTS; e.g., Schiissler et 
al. (1996)]. The fact that the fields of T Tauri stars are not ax- 
isymmetric may have important implications for the accretion 
process. Our work suggests that dipolar fields tend to divert a 
significant fraction of disc matter into the wind. Channelling 
the disc material along magnetospheric stellar field lines be- 
comes more efficient if the magnetospheric accretion process 
happens in an episodic fashion. It is conceivable that a nonax- 
isymmetric stellar field has a similar effect, but it would then 
produce time variability on a time scale similar to the stellar 
rotation period. Variability on time scales similar to the stellar 
orbital period has indeed been observed; see Johns-Krull et al. 
(1999). 

We find that the sum of magnetic, accretion and viscous 
torques acting on the star, is positive at most radii between the 
star and the disc if the stellar surface field strength exceeds 



a certain value that is somewhere between 200 G and 1 kG. 
This positive torque means a stellar spin-up, in agreement with 
Romano va et al. (2002) for stellar rotation periods of around 
9 to 10 days. The accretion and viscous torques are almost al- 
ways much smaller (in amplitude) than the magnetic torque. 
The small accretion torque in our simulations might be due to 
the relatively small accretion rates and accretion flow veloci- 
ties, both of which are not imposed. These, in turn, might be 
due to the fixed inner disc radius. 

The mean-field disc dynamo is responsible for the struc- 
ture in the disc wind, with coexisting pressure-driving of the 
outer disc wind and magneto-centrifugal launching and accel- 
eration of the inner disc wind. The stellar wind is always mostly 
pressure-driven. 

In contrast to the simulations of Hirose et al. (1997) and 
Goodson et al. (1997), we do not find any stellar (fast) jet 
that is driven by magneto-centrifugal processes. However, our 
magneto-centrifugally driven fast inner disc wind is close to the 
star and might collimate at larger heights to form the observed 
protostellar jet (cf. Fendt et al. 1995). 

Another important result that has emerged from the present 
investigations is an anti-alignment of the disc magnetic field 
relative to the stellar dipole (no X-point). In our models, the 
relative orientation of the disc field and the stellar field is no 
longer a free input parameter, but a result of the simulations. 
The only obvious way to prevent anti-alignment of the disc 
magnetic field, is to have magnetic fields in the disc and in 
the protostar that are entirely due to the accretion of an am- 
bient large scale field. This may indeed be quite plausible for 
many protostellar discs in a large fraction of star forming re- 
gions. Another argument in favour of this possibility is the fact 
that strong collimation of protostellar outflows into well pro- 
nounced jets has so far only been found in the presence of a 
large scale field aligned with the rotation axis of the star-disc 
system (e.g., Ouyed et al. 1997), and not for dynamo-generated 
disc magnetic fields (see Paper I). 

An important limitation of the present work is the re- 
striction to axisymmetric magnetospheric accretion. Given that 
fully three-dimensional simulations have now become feasible 
(e.g., Hawley 2000), this would certainly be an important con- 
straint to be relaxed in future simulations. 

Appendix A: Effective stellar rotation rate 

In this appendix we give a precise definition of the effective 
stellar rotation rate, O*, as it was used in Eq. i2l\ . 

Integrating Eq. ^ over the stellar volume yields an evolu- 
tion equation for the total stellar angular momentum, L„, 

^* ^acc ~T" 7/' ma g -|- T v [ sc , (A.l) 

where L* = J gl dV is the integral of the angular momentum 
density, gl, over the stellar volume. Obviously, can grow 
from mass accretion alone without spinning up the star. We 
therefore need to look at the evolution of the specific stellar 
angular momentum, I* = L*/M*, where M* = J gdV is the 
mass of the star. Using the product rule, we can write the left 
hand side of Eq. JA. ll > as 

L» = /'*M» + /*M», (A.2) 
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where the accretion rate M* can be expressed in terms of the 
mass flux density, gu po \, via M* = — § gu po \ ■ dS, where 
the integral is taken over the stellar surface (dS is the out- 
ward directed surface element). Further, Z* = L*/M* = 
J gldV/ J gdV can be expressed as Z* = OT*f2*, where 
is defined as a weighed average, 

w 2 = (vj 2 ) = <j> gu va \uj 2 ■ dS I <j> gu poi • dS, (A.3) 

where the integrals are taken over the stellar surface. This de- 
fines O*, which replaces the intuitive definition of fi* as the 
stellar surface rotation rate (the latter, however, is still suffi- 
ciently accurate a definition for all practical purposes). 

With this definition of tu* (and f2*), we obtain an evolution 
equation for the specific stellar angular momentum, Z*, 

A^Z* ^acc Z^A/* ^inag ~t~ -^visc; (A. 4) 

which can be written as 

M*l* — ^acc,nct ^mag ^visc (A. 5) 

with T acc nGt given in Eq. (121 1 . This shows that the sign of the 
right hand side of Eq. iA.5\ determines whether the star spins 
up or down. 
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